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Abstract Based on the concepts of infinitesimal and super real time-space composed 
of monads, the equations describing the particle flow are proposed. After the initial 
conditions are given, the state of the particle system at any time can be obtained by 
solving these equations. The physical properties of the granular matter, such as the 
moving velocity, the rotation angle, the angular velocity of rotation, the temperature, 
the pressure, the density, the mass, etc. at any time and space position are obtained. It 
is also to obtain that the translational motion is not dependent on rotation motion of 
the particles and can be separately analyzed. For a stationary particle system, the 
conditions that must be met by the force and force moment on the particles in the 
system are given. On this basis, the properties of the granular matter at rest can be 
further explored for. 
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1 Introduction 


Granular matter is common in nature, daily life and production technology, such 
as sand, soil, ice floes, snow in nature, food, sugar, salt in daily life, coal, ore, 
building materials in production technology. Some medicines and chemicals are also 
granular matter. There are also many discrete material systems such as bulk cargo 
transport, earth plate movement and the flow of vehicles on the road. They are also 
often treated as particle systems. It can be said that granular matter is the most 
abundant on the earth, one of the most familiar types of substances. 

The annual production of cereals and other granular matter in the world is in the 


order of tens of billions of tons. The production, transportation, processing and 
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storage of these materials consume about 10% of the energy of the earth each year. A 
deep understanding of granular matter will bring great benefits to the development of 
global industry and economy. At the same time, in order to prevent natural disasters 
such as debris flow, avalanche, ice floe, landslide, desertification and earthquake, it is 
also necessary to understand the movement law of granular matter. 

For centuries, many famous scientists have studied the movement of granular 
matter. Coulomb [1] first proposed the relationship between the slope angle of the 
sand and the friction coefficient. Faraday (1831)[2] found that the vibration causes the 
particles to form convection and accumulate , Reynolds (1885)[3] pointed out that if 
the particles are piled up in an elastic bag, any additional action will increase the 
volume of the particle system, which is called the Reynolds principle of pressure 
expansion. British scientist Roberts (1884) [4] first noticed the "granary effect" of 
granular matter. 

Although, in order to applications, the industrial, engineering, and technical 
communities have studied the production, processing, and transportation of granular 
matter for a long time, the rise of physicists' interest in research of granular material 
has been nearly two or three decades. People have carried out experiments and 
simulations on various aspects of granular matter, and obtained many meaningful 
results. But their understanding of the laws of motion is still superficial. The basic 
theory of describing granular matter has not yet been established. Some problem of 
the most basic aspects of granular matter is still plaguing people. Regarding the 
understanding of granular matter, the famous theoretical physicist de Gennes [5] 
thinks that almost everything in this field remains to be understood. The famous 
physicist Kadanoff [6]said that ordinary fluid equations cannot be used to describe 
granular matter, and its rich and peculiar behavior is not well understood. 

Granular material is a complex system. The study of granular matter not only has 
an important application background, but also a deep understanding of its state and 
motion law will promote the new development of physics. The work in this paper is 
based on the concepts of infinitesimal and super real time-space. And in the work a 


method to describe the movement of granular matter is proposed and the equations 
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describing the law of the particle flow are obtained. 


2 Description of physical phenomena 
and its mathematical basis 
In order to describe the physical phenomenon in which the properties of physical 


objects change with time and space, the current physical practice is to divide the 
objects into small units in time and space Av,At , and the properties of the objects in 
the small range are considered to be uniform. 

Therefore, knowing the value of the physical quantity on every range Av,A? , the 


description completely of the physical quantity of the entire object with time and 


space is obtained. That is just the description of the physical properties of the object. 
Here, Ax,At is the "sufficiently small" unit, the meaning of "sufficiently small" is 


small enough to reflect the changes of physical quantity with time and space. 
The mathematical abstraction of such a description method can be summarized as 
that the object is divided into many small units, each unit has a spatial scale Ax ; the 


time size of the physical problem is also divided into many small units, and the scale 


isAt. Then let Ax —>0,At—>0 i.e. tend to time-space point (x,t) to define the 
value of each function f representing physical quantity at the point (x,t), and 


further to define the various order time and spatial partial derivatives of each physical 


quantity function f in the physical problem, such as: 


Opa Oi Of 
ôt ôx O°t  A’x 


etc. Their definitions are well known, for example 


F un EAD- SO 


Ot ^0 At 
E e (1) 
Ox x0 Ax 


On the basis of (1), the high order partial derivatives of the function can be also 
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defined. And the law of the motion properties of the object with time and space is 
described by the differential equations containing these derivatives. 
Mathematically, definition (1) is based on the concept of limits (called 


-e framework). In order to understand the mathematics meaning of (1), we must 


carefully examine the meaning of Ax —> 0, At + 0. It is very clear to the physics 


researchers that Ax > 0,At—0_ in (1), in physics, does not actually mean that it 
tends to be absolute zero, but rather that the physical properties of the object are 
uniform in a sufficiently small range Ax, At (or the variation in this small range is 
negligible). 

That is to say, when the physical properties of the object are uniform in a 
sufficiently small range Ax, At , the mathematical analysis used to describe the motion 


is based on (1), i.e. based on the limit concept under the -e framework. And 
when the motion law of the object can be expressed by the mathematical equations 


based on (1), the motion properties of the object must be uniform over a sufficiently 
small range Ax,At (or its variation within the range Ax, A? is negligible). 
However, in the physical world, not all motions are that in a sufficiently small 


range Ax,At the physical properties of motion are uniform. Some people may think 


that if the physical properties of the motion of the object in a range Ax,At are not 


uniform, the segmentation can be made smaller, and a smaller one can always be 
obtained, so that the physical property of the motion in this smaller range is uniform. 
In fact, this opinion is not appropriate. 

Because such a statement ignores a key fact that when describing the physical 


properties of an object's motion as a function of time and space, people divide time 
and space into many small units Ax,At, and the size of the small unit Ax, At is 
corresponding with the specific physical problem studied. That is to say, it cannot be 


subjectively and randomly segmented by people. For example, the small unit scale in 


the flow around a cylinder in a trough is clearly quite different from the small unit 
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scale in atmospheric physics. If people leave the specific physical problem and 
randomly select the scale of the small unit, there will be difficulties in physical 
description, such as excessive calculation. 

Since in the specific physical problem, the scale of the small unit can be not taken 
subjectively small, but the small unit scale has objectivity, then it is possible that, in 


some physical problems, the physical properties of the object motion in the small unit 
range Ax,At are not uniform. As it happens, turbulence, particle flows are such 
physical problems. The discontinuous medium of the particle assembly obviously has 
a non-uniform physical property in the small unit range Ax, At . 


Scientific description of such physical problems is unreasonable in mathematics 
based still on (1). It is necessary to propose a mathematical concept that satisfies the 


variation of the values of the functions representing physical quantities within the 
small unit Av, At. This mathematics is a non-standard analysis proposed by American 


mathematician A. Robinson[7] in the 1960s. 

In non-standard analysis, there are non-standard numbers in addition to the 
standard numbers, where infinitesimal € is an important non-standard number. 
Positive € is a number greater than zero but less than any positive real number. All 


real numbers are called standard numbers. Infinity ZL is also a non-standard 
number, a+ €,aé,€€,LL ..- all are non-standard numbers (where @ is any real 


number). 


In the standard analysis, any time and space point can be expressed as 
(xX X%3;t ) in a certain coordinate system, where ,,*,,X;,¢ all are real numbers. 


These time and space points are geometric points, which are absolute geometric 
points with zero volume, usually called real points. For convenience of explanation, 
take one-dimensional space as an example. As is well known, the real number 
one-to-one corresponds to the point on the real axis of the one-dimensional one, i.e. 
any point A on the real axis corresponds to a real number @ , whereas any real 


number a must correspond to a point on the real axis (Fig.l). However, 
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non-standard numbers cannot be matched by points on the one-dimensional real axis, 


and points on the real axis cannot be used to represent non-standard numbers. 


A 
—e y 
Figure 1 Real axis 

For expressing the non-standard number, the real point A corresponding to real 
number @ should be expanded into a monad, which is a collection of numbers. Any 
one of monads contains only one real number @, and the remaining numbers in the 
monad are non-standard numbers. This monad is called @ -monad. The scale of the 
monad is infinitely small (Fig. 2). 

Thus, the real number point is expanded to a monad and the real number axis is 
expanded to a super real number axis. The real number axis is composed of real points, 
which are absolute geometric points whose scale is zero. The ones which the super 
real axis is made from are monads. The monad is not an absolute geometric point. The 
scale of monad is infinitely small € , which is greater than zero. Since the scale of 
an absolute geometric point is zero, the absolute geometric point has logically no 
interior. If you would like to take a physical small area as an absolute geometric point 
in mathematics, it is only reasonable that the physical quantity in this physical small 
area is uniform. However, if a physical small time-space region is mathematically 
abstracted as a monad, the scale (volume) of the monad is not zero but greater than 
zero, it is logically reasonable to say that the monad has its interior. Then the physical 
quantities are allowed to vary and be not uniform in this small region, i.e. in the 


monad mathematically. 


onad 


Figure 2. Monad (A stands for real number a ) 
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Through the above analysis, we know that when the physical quantities in the 
small range Ax, At are not uniform and vary, our description of the physical 


phenomena cannot be based on the standard analysis of (1), but based on non-standard 


analytical mathematics. 


C aa o oS 
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Figure 3 Hierarchy 


It is easy to know that the monad is a lower level time-space. And Figure 3 
shows three levels. The B level is the 0 -monad in the 4 level, the B by 
comparison to A is the lower level; the C level is the 0-monad in theB level, 


which is one level lower than the B level. For the convenience of description, a 
coordinate system (x',t') is established in the monad, x’ parallel to x (same 
direction) (Fig.2). 
So for both(x,t) and (x’,t') levels, there are 
dx=eé,dt=e, dx'=6,dt'=e, (2) 
But from a hierarchical perspective, the true lengths of (dx',dt')in the level (x,t) 


are(here ao are the dimensionless infinitely small): 


o2 o2 o2 o2 3 
, , 
E dX =E € Edt =€: £, (3) 


Take a point x in the real number axis. The numbers inside the x -monad can be 


represented by coordinates x’, called a internal point of the ¥-monad. And use a 
function f(x,x’) to represent the function value on the inner point x’ of 


the x -monad. Can also be defined 


of (ater) B f(xte,tx, t')-f (x,t,x',0') (4) 


Ox g 


x 


Of (x, t, x',t') f(x tt+e,,x,t)-f (xt, at | 

SS eee (5) 
Ot E, 

For the four-dimensional space, there is similar discussion, i.e. the real point 


(xi X3, X3) is expanded into the monad (xi X2 X3) , containing only one real 
number point (%,5%)5X3.t) . (X1 X2:X3;%,X1,X3X3,t') denotes the internal point 
(x{,x5,x5,t') of the monad (x,,x,,x;,¢). The function value of the internal point 
(x{,x5,x5,t') of the monad (%,,x,%;,f) is represented by a function 


ror ngs 
J (XX3 Xb, X1, X3 X3 ), and so on. 


In the standard analysis, the variables in the equations describing the motion law 
of the physical object based on (1) are the functions on the absolute point. In the 
non-standard analysis, the description of the physical motion law also requires the 
equations. In this case, the variables in the equations based on (4) (5) are often related 


to the whole monad, such as some average value. 


3 Description of granular motion 
(model of granular motion) 


The time and space in which the granular matter is located is divided into small 
units Ax, At. Particle motion is a physical phenomenon in which the motion 
properties (physical quantities) are variable and non-uniform in the small unit 


range Ax, At. For example, the two small spheres are close to the point of contact, 
there is arbitrarily small Ax,At around the contact point . And the inside of Ax, At 


is non-uniform, there must be a space of density 9 #9 and space of density P=0 (i.e., 


no matter) (Fig. 4). Therefore, in mathematical abstraction, the movement of granular 
matter occurs in a super real time-space composed of the monads. That is, the particle 


matter flows out from a monad and flows into an adjacent monad. 
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Figure 4 Near the contact point of two ball 


Regarding the translational motion of particles (or particle parts), they can be 
classified. Let there be a total of k particles (or particle parts) in the monad 
(Gar 9% ) . They are marked as 1,2,---k . And their movements can be classified as 
three types (Fig.5): 

A. Movement purely inside the monad, such as particle "1"; 

B. The particles are on both sides of the interface between the adjacent monads. 
Such as particle "2", "3", "4". Particle "2" is flowing out of the monad, "3" 
flowing into the monad, and "4" can be regarded as "2", "3" in combination of two 


situations, both flowing into and out of the monad; 
C. The movement of particles (or parts of particles) at some time, such as t =t; , 
occurs with inflow (or outflow) of the monad. For example, particle “5” is outside the 


monad before the time t; , and becomes “3” type after the time ti . About the 
movement of particle "6", beforethe t; , there are some materials of the particle 


in the monad, and at the time t; particle "6" completely leaves the monad. In fact, 


the (C) type of motion is the movement being adjacent to starting and ending of the (B) 
type motion. 

The parts of the same particle in the same space-monad have the same physical 
properties; while the physical properties of different particle in the same space-monad 


are different, which reflects the discontinuity of particle system motion. However, the 
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physical properties of the same particle divided into two adjacent monads have the 
infinitesimal differences, which is consistent with the continuous motion of the 


material and is the continuous side of the particle flow. 


Figure 5 Classification of particle translational motions 


On the particle flow, a certain limit is imposed on the particle characteristic size in 
the system. If the particle size € (monad scale), the flow of a large number of such 
particles becomes a continuous medium motion, which can be solved by using of the 
hydrodynamic equations. 

Now the particle flow we are going to study is not the above case. It is now to 


consider the system composed of a large number of particles, while the particle 
characteristic size ~€ or greater than ¢ (monad scale). The movement of the 


aggregate system composed of a large number of such particles is called the particle 
flow or the movement of granular matter. What we want to describe is the movement 
of this situation. 

For the sake of simplicity, it is assumed that the particle matter is an elastic 
substance, and the influence of a fluid such as air on the movement of the particles is 
ignored. 

Here are a few useful concepts. 


(1) The inner point of the monad 
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The super-real space is composed of the monads, and the monad is composed of 
the lower-level monads. The lower-level monad can be considered as a geometric 
point in a sense, that is, the physical properties inside the lower level monad are 


uniform. This low-level monad is called the inner point of a high-level monad. 


The scales of the monad are infinitely small quantities €,,€,,&3,€,, here &; is 


infinitely small in i-direction of space, €&, is the time infinitely small. The scale of 


. oO 2 oO 2 oO 2 oO 2 . 
the lower level monad is E, 61 ,€) 62,8, 63,E,Er * The volume of space monad is 


2 
oO oO oO 
€,€,€; , and the volume of lower-level space monad is (a E2 3 CEE o 


Therefore, the number of inner-points in the space monad is 


-2 
EE E oO oO oO 
N, =——" =| @1 £28: 
eo (6) 
EEE | E1 E2 E3 
Similarly the number of inner-points in the time-space monad is 
EEEE bo oo y” 
N! = Laca =| £1 £283 E1 
(7) 


oO oO oO oO 2 
EEEE | 6162 €3 Et 
o o $ $ ä + oO o . . 
Among them, ¢;,¢, are the dimensionless infinitely small, 7, 7 dimensionless 


infinite. 


Take 
Liéi =1, Té,=1 (8) 
(2) Particle and the element of a particle 


System is composed of a large number of particles, its motion is the object of our 


research; the particles move between the monads. And we call the particle material at 


an inner point of the space monad as the element of a particle. Since the inner point of 


the monad is a lower level monad, we believe that the particle material is evenly filled 
in the inner point of the space monad. That is, the density of the particle material has 


only one value for one inner point, i.e. the density at the inner point is uniform. And 
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so are other physical properties. 
In space, the element of a particle is always at an inner point of the space monad, 
moving from one inner point to another. 


(3) The mass inner point and no mass inner point 


Ata time-space (£1, £3, €3, & ) , the total number of inner-points is Nj} . Any one 


of these Ni inner points which is not occupied by the particle matter is called as 
no-mass inner point. Otherwise, it is called as the mass inner point. 
Between t —>t+e€, , if an element of particle moves inside a monad (including 


stopping at a space inner point), the number of mass inner points contributed by the 


element of particle is: 


goa oe (9) 
EF; 


If the time for which the element of particle stays in the monad 


is(t, t=0 3914, = t') , the number of mass inner points contributed to the monad 


by the element of particle is: 


5. 2 


UE M 
a (10) 


oO 


Ee; 
An element of particle flows out of the monad, and there is a corresponding 
number of mass inner points in the monad flowing out of the monad; when an element 
of particle flows into the monad, there is a corresponding number of mass inner points 


into the monad also. And let an element of a particle flow out of(into) a monad at 
(t,t') , the corresponding number of mass inner points flowing out of(into) the monad 


is 


R (11) 
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Let there be AN elements of particle flowing out of the monad (x+¢) (Fig. 6). 


They flow out at interval of the(t,t' =0)> (t, =a’) . 


X X+E 


Figure 6 Elements of particle flowing between monads 


Then at interval of the (t, t' =0) > (t,t = a’) the total number of the mass inner 


points flowing out from the monad (x + é) is 
AN' = f ——dn (12) 


In case of equal sections and equal velocity of elements of particle, you 
have dn=adt' ,here a is a constant. 


Then 
AN' = =< AN ( ) 


Generally speaking, there is: | AN’=J(AN) (14) 


(4)The average over the monad and the average over the substance of particle 


Taking the physical quantity f (x1, X3, X3,t, X1, X3, X3,t'), its average over monad is 


~ TS pee a ae 
f =f t) =7 [ax GS | ax, zjar (15) 
3 0 0 


Let B be a collection of elements of particles (a collection of particle substances), 


and B's contribution to the average over the monad is 
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~ 1 
f tee 2 = 


aa | faxidx,dx;dt' (16) 
172-3 SB 


There is also average value as follows: 


~B ~B 


1 1 d (4 i 
f =f (anh B= Fey TREA (17) 
B 


V(B) is time-space volume of B . This is the average value of f over B, 


called as the average value of f over the collection of related substances in 


time-space. 


It is available from (16), (17) that 


NGF (xo xX,t,B)= N'(B)F (x,.25.%55t,B) (18) 
Where N’ (B ) is the number of mass inner points by B . 


(5) Strain, stress and elastic force in the particle 
The space in which the particle system moves is composed of the space monads, 


and the movement of the particles is between the monads. The part of the particle / 


in the space monad(x,,x,,x;) is marked as(x,,%,,%;,/) (or/(x,,x),x;)), and its 
average velocity over the substances is U, (eshte al ) . The velocities of all elements 


of particle / in the same monad are U, (tadesmitl) f 


According to the usual practice of elastic mechanics, the strain and stress of the 
particle / can be defined as follows. 


Displacement change 


2 


Ar; Seite. +E Xbt L)-U, (xx, Xbt) |e, i (19) 
Strain 

o2 
e(tr)= 2 =[U, (x +6). jx stat sl) -U, (2.2% stat st) PO (20) 


L t 
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7; (t,') =. Shey E #6 otatit d) U(x pot) E 
o2 
+U, (x, $6 kata d)—U, (ox atl) JE (21) 
Consider pre-strain 
e, (t,t') = ef (4,0) +e; (4,0) (22) 
rir)=ri (G0) +7) (60) (23) 
Similarly there is 
0 =e, +e, +e; (24) 
Then the stress is 
o, =40+2ue, (25) 
Tj = HY; (26) 


Here À, 4 are Lame constants. 


The elastic force is equal to the stress multiplied by the corresponding area. 
Due to the discontinuity, this corresponding area must be defined and sought for. 


The area is given below. For example, the particle? moves through plane A-A. The 
number of elements of particle passing through the plane A-A within (t >t’ +dt') 


is An, the average number of elements of particle passing through A-A in unit time 


An 
is ,2. And the time interval, in which one element of particle passes A-A , is 
dt' €; 

o2 
E Ei 


A Hence the number of elements of particle passing through A-A during 
U, (tt 


that time interval is 


2 


E, a An 
= (27) 
Tl) yk 


These elements of particle are just the elements of particle in the /-cross section on 


the A-A plane. This cross section consists of these elements of particle (a layer of 
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elements of particle). Obviously the area of this layer of elements of the particle is 


o2 


. o2 a o2 : 
s(t +d)= =a as E — ar as (28) 
(Bt) we i( 2 P 


Note that the values of the physical quantities (density, average velocity) on two 
sides of A-A plane may be different, so there are two areas on two sides of A-A plane . 


We take the smaller of the two as the area in the calculation of the elastic force. 


4 One kinematic condition 
Let U;(%,.%25%3,,1,%5,%5,t') be the velocity of movement of the granular 
matter (which consists of a large number of elements of particle, and the element 
of particle is indicated as (;,x5,%;,t’)). Its average over monad is 


1 
“EET 


U, (xX X3,X33t) ff ff U dxidx,dx;dtť 


Put case that there is only one element of particle moving in the monad. Then 


have 


E 


1 POP l „E; 
z J Ui (eae grle 


7 1 d f 
== Jue jdt'e,6,€, = — 


1 
o o o o =l 
L (a E2 sa ) 
(29) 


d 
Among them, f dx; represents the integral along the motion trajectory of the 


oO 2 oO 2 
element of particle, and the relationship is also used here: [7 (t')dt'e, =dx'e, 


Then there is, U, = rere —x!,) (30) 


Here *;..%;, are the lower-level coordinates of the end point and beginning point of 
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the motion trajectory of the element of particle respectively. 
(30) is the contribution of only one element motion to the average over the monad. 
There is now an assembly of elements of particle B in the monad, which contributes 


to the average over the monad by 


FT 1 i ' ' 
OR ERA E E(x) (31) 
NT 2° 


t 
Here, the meaning of 2 is to sum over all the elements of particle in the B . 


By (18), there is 


NOU (x%,x3,%,t, B) = N'(B)O (Get, 8) (32) 


N' (B ) is the number of mass inner points within the assembly B . Then there is 


l i ' ' 
U (2% %ytB)=— a > (4, -%,) (33) 


o N, € B 


1 surface Į surface 


Figure 7 Definition of surface i°,i 

For the space monad ETA , Figure 7 shows the x;- direction. Then the 
surface x; is taken as the surface ; and the surface x,+€, as the surface į". For 
a monad (t Xe ) , the coordinates of the start or end points of one element motion 


trajectory on the surface į are all zero; the coordinates of the start point or the end 
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point on the surface i* are infinity L, . The start point of one element of particle is 
located on the surface j*, it is considered to flow into the monad through the surface 
it. The element end point is located on the surface j*, it is considered that the 


element is flowing out through the surface j* , and the rest direction is analogous. 


The right summation in (33) can be written as 


Die =x!) = > x = > Xii + N” L, = NUL, (34) 


x} .(B)fin. x} ,(B) fin. 


Then (33) can be written as (by using a Ta zi) 


N'(B) a 1 E, i 
Clie AN AB E a x= xi, +N“ - N" 
N! i ( 1929-3 ) pa NT pe Z | ae NT ( i i ) 
Et Et 
(35) 
Here , on - 2 are the sums of the terms having finite coordinates of the end 


point and beginning point of elements of particle respectively. And N 3 N i are the 
numbers of the elements of particle flowing out and in the monad through the surface 


i* respectively. 
It is easily to know that only (35) is independent kinematic condition. 


Let m, , related to density, be the mass of one element of particle, ™,,m, be 


the mass of the B at end time and the starting time respectively, x;~,x;), be the 


coordinates of the end point and the beginning point of the centroid of the B, then 


there is 
o2 
N' B i 1 e e E; 1 ou. in 
y ly EEATT] = = N Tm, (mxi -mxi )+ A NT (N Ni. ) (36) 
t t 
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5 Translational movement (A) particle interaction 


First there is the discussion of the motion of the "2" type in Fig. 5. 


The particle / leaves the monad (x, a and enters the 
monad (x, —£,,X;,%;) (Fig.8) . Let the velocityU, (x; (a) say (a) sm (q).t,',1) be 
given. Among them, (q ) represents six adjacent monads or the six interfaces of the monad 
(eh), here q = (1,2,---6) or q =(, 7k ,i°,j*.k*). And 
(x5 Xps Xps X X’, x) will be abbreviated to(x,x'). And following the steps below , the 


translational motion of / (xax ae | in the monad can be solved. 


(x, -,,x,,%,) (atn) 


Figure 8 Particle / moving in the type of "2" 


(1) By the velocity U, (x; (q4)»x, (4), x, (4),t,t',1) in the adjacent monads, 
from the relationships (25), (26) the elastic deformation stress in / , 
ie., (x,t,t',1)  Tp(x,t,ť',l) , are obtained. Because s;(x,t,t',/) is known thereby 
further the elastic force F,(x,t,1',/) is obtained. 
(2) The total force on / (xx Ea exerted by the other particles (or particle parts) 
in the monad a) is f, LERET A , then there is 
LiFe (ox tott) =0 (37) 
(3) According to the conservation of momentum, there is a relationship 
|F, (BE, (4,2 ).%.6t52) dee 
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= U, (xp xxt t + dt’) m(x,,x,.%,.60 + dt’) 
-U, EEEF IREE 7) 


+8m™ (e = i +dt',1)U, (oxx) 


-Y Sm” Cee ad > t'+dt',1)U,(x,(q),%, (q).%;, (q).t,t'.1) 


q=1 


Ignoring high-order small quantities, there is: 


|r, (x,.%;.Xpstst lL) +f, (oxp ms6t5) dee 
=| U, (xx; Xt, t +d’) —U, (x,.%).%, 6.001) |m(x,.%).% 60 +dt',1) (38) 


Among them 


m(x,.x).X,t./ +dt’,1) = m(x x xpt t eon (aa 1’ +dr',1) 


Sm" (x,.%) Xt — 1’ +dr',1) (39) 


q=l q=l 
By (38) can obtained 
F, (iat (ttad) ~m~ E? (40) 


(4)Now give the method to get at f(x ,%.605/) , 

et Fe ee (41) 
When, for a certain particle 7! , F,(x,x,.%,60,/)=0 , 

bet f(a tt Tamla ett) (42) 


Here the [a,] is acceleration dimension. 
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k 


And Dr ae) =i), (oita) (43) 


l=1 


Then ARET will be allocated by size of Tord) „i.e. let 


k 


K Cree eee Ža, =a 


l=1 


and take eaae ced . (44) 
a 
k 
Obviously Ila )= ARER (45) 
fal 


Here d, is a positive constant. 


Then take the amendment to f” 


AER ETA =f" (at-a, (oaa) (46) 


Obviously there is YL. (riety ttt!) =0 (47) 
m 

That is, AE is corrected to AOTER by 
A EEE . Such Tete.) obtained can be taken as the total 
force on / (es x TA exerted by other particles (or particle parts) within the monad . 
(5)The velocity U, CEER at ) will be obtained when substituting 
F.(%,%;%pste's1) ino G8): 

And in the formula (38), m(x;x;,x,t,t'+dt',1) is required, which will be 


discussed in the next section. 


(6) Repeat the above steps to find the velocity of EA pt and the internal 


forces (the total force on / (xax Be | exerted by the other particles in the monad) at 


any time. 
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About interior motion purely ("1" type in Fig.5) 


With the interior motion purely of the particle /' (x,, X; Xystat ) , known the initial 
velocity U, (oaia) ’ let 


ia mi ac el 1) a, 
Then to correct f? (aaa (7=1,2,---I'---k) and get at each 
Í, (etaa) , included f, aa et 1) . So the acceleration of the 
internal motion a, (aa Xpt, of ‘) can be found, and so the next moment velocity 


U, (atant a) . In order to proceed to the next step, let 
f (totti adl) =f, CERTA 


Repeatedly you can get at the velocity U, (xx a) and internal force on /' 


at any time. 


About motion in the "5" type (Fig.5) 


The particle "5" completely outside the studied monad Cae enters the 
monad at some time between (t '=057=T ) . Let there be not “5” in the monad 


(wax oe before the timet/=t; , while t'=t; particle “5” starts to enter the monad 


CEEA from the(x, -EXX ) monad. Because, at the timet’=¢,, the mass of 
particle "5" is very small in the monad (1.855%) , it can be set 
Hat etot 1,5) =U, (ee, ott =S Thats, 

AEA EE 
So let fees) = e l] 
Then by the amendment to f” , you can get at f, (Cad). by (38) the 
velocity lutan Ue 5) will be obtained. And further 
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F, (x, Xis Xy»t, t +dt', 5) can be obtained. Then you can follow the previous method. 


6 Translational movement (B) Formulation of Equations 
Assuming that the particle / flows out of the monad (x, a) and enters the 
monad (a = Eee) (Fig.9). Between t' > t'+dt' , the part of the flow out of the 
monad (x,, x ee, (flow into the monad( «x, =E X; X; ) ) is denoted as A, and the 
mass of the outflow is mp” e >t'+ dt',1) , and likewise, the part of the 
flow into the monad (x, Sak on (from the monad (x, Xe )) is denoted as A’ ,the 


mass of the inflow into the monad (x, EN a) is 


ôm; (x, ~6,,X,,%,60 >t +dt',1), 


x.-6&, x 


i i i 
Figure 9 Movement of the particle 7 between two monads 


The unknown quantities required are the quantities at the next time ¢'+dt’. For 


the /— movement between two monads the unknown quantities are as below. 
On I(x =E, X aw) there are pressure, density, temperature, velocity, 1.e. 
P(x, -Ep X Xp bt +dt',1), p(x, —€,,X;,%,5t,t' + dtl), 
T (x, -EX Xpt, t +dt',1) ,U, (x, — &,,X;X,5t,t' + dt',1) 
On / (ae E there are also pressure, density, temperature, velocity, i.e. 
P(X;.X).Xpstt'+d0',1) f p(x xX; xXyt,t'+dt',1) : 
T sarab aI) U; (xxx t +dt',1) 


Moreover, N’ H (4') , N an (A) . So there are total fourteen quantities unknown. 
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The equations for solving these unknowns are given below. 

The equations for / (x, SEX j Xy ) are: the energy equation, the momentum equations, the 
equation of state and the pressure formula at t'+dt'. 

The equations for / (xs a) are also: the energy equation, the momentum equations, the 


equation of state and the pressure formula at t'+ dt’. 
Moreover there are the energy equation of 4 and the mass conservation equation for 4. 
There are 14 equations above, so the problem can be solved. 


Now list these 14 equations below. 

The energy equation on / (xx =A 
m(x,,X,,Xpst +dt',])u( xxx 60+ dt l)—m(x,.x Xt u(x xxt) 
-őm(x;, x, xpt, t' > t'+dt'1)u(x,,x,.x, 0.01) =()" EEN > 1 +dt',1) 


E m\ tbt +dt' 1) m(x,x,.%,6051) 
DP Megat) EEIT] E T (48) 


The momentum equation on / (ax A 
| F, E E E 
=[U, (xxx tt tat l)—U, (x,.x).2, 0000) |m(x,x,, X136, +dt',1) (49) 


The state equation on Iltut at) at t'+dt' 


1 1 1 mat na ' 1 
EEEE, = DATO) [Leal T(x.x)a 60 +dt',1) T,) EEIE EE, +dt',!)| 


(50) 


The pressure formula on /(x,,x,,X,) at ¢/+dr' 
P(x) ,A stl bat’) =—2 (040m +0)( Ht teste + dll) (51) 
The energy equation on / (x, = Cy X; x) 
m(x% -EX Xpt t tt’) u(x, -Exit t tat’) —m(x, -Exx t 1)u(x, Ekt, t) 


-ôm(x, -Ep X Xbt ->t +dt',l)u(x, Ept] =0" (x, -Ep re b, t'+dt',1) 
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E 7 m(x,—€,,%),% A] mx, —€,,%),% 6,01) 
p(x EX pXatt'l) P(x; -Enx Xt,t'+dt',1) MEETER ne 


The momentum equation on/ (x, 65% 3%; ) 


2 


| F, (x, =~ 6,,X;,%p stl’ l)+ S, (x = 65%; %patt' 1) | 
=|U, (x, -8X Xpt, t +dt',1)-U, (x, — 6 5; Xbt,1) |m(x, -8X Xp t,t'+dt',1) (53) 


The state equation on Hone x 2; } at t'+dt' 


1 1 1 7 d | 
me eer sy) =T) [i+a(T(% -8th +dt',l)-T,)—Kp(x,-&,.x,.%stst +dt',1)| 


(54) 


The pressure formula on I(x, = Skai] at t'+dt’ 


p(x, = €,,X),X, t,t + dt’) =-3(6% +05, +033)(x, —£,,X,,%,.,U +dt’,1) (55) 
The energy equation on A : 
[u(t,0' +t’, A')-u(t,t', A) |m( A) 


=-p(tfť', ) l 
p(t.t'+dt',A’) p(t,t,A) 


m( A)+0" (t,t! >t +dt', A) (56) 


The mass conservation equation: 


p(x -EX .%p ott tdt I) NË (x, -85x Xbt ot +dt'1)V, 
=p it prab t NE okt >t NT, (57) 
Where u is the internal energy per unit mass, which can be obtained from the equation of 


state in thermodynamics as a function of the temperature and density. Q” is the input heat. V, is 


the volume of internal points of a space monad. The energy equation is just the first law of 


thermodynamics. 


The cross section s; (x, t,t'+dtť',l),s (x=e,t,0'+ dt',1) can be obtained from (28) , and 


will be used for the next step calculation. 
In addition, kinematic conditions are required for the position of the centroid. 
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Applying kinematic conditions to the outflow part through the face i of the monad 


(x Xis Xy ) (coincident with the face i” of the monad (x, — Ej,X Xy ) ) A, you can get: 


N' A f i Ic out 
We sit) AS (ANN (4) (58) 


And applying kinematic conditions to the part A’, there is 


1 ' o2 
MAy (x -expt t tal) = fa se )-«,} No) (59) 
g NTE 


From (14), we can see the relationship between N '(A) and N” (A) (also N '(A') and 


By (58). (59) xi5,(A) and xs (4) can be obtained. The coordinates in the direction 
X, X, of the centroid A, Acan be processed as follows. Let the components in the direction 


X;,X, of the centroid-coordinate of the cross section of l at aiea] on the interface 


te 


“i a respectively (the subscript g indicates the 


between the correspondent monads be x 
interface between the monads), and obtained from the pre-determined initial conditions. It is 


possible to make the components in the direction x,,x, of the centroid-coordinate of A be 


ii (Pak at the time (t „t '), so that the components in the direction x,,x, of the centroid- 


coordinate of A’ at the time (t „t'+dt ") can be easily obtained. With the centroid- coordinates 
of A ad A’, you can calculate the centroid position of 


I(x; x; Xpt,t' + dt’ )and/ (x, ~6,,X;.%, 6,0 +d’). 


7 About the rotation of particles 
1, Rotational angle of particle and setting of some coordinate systems 
If there is the deformation inside the particle, there is the relative displacement 


between the points inside the particle and the particle state is not equilibrium. In order 
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to describe the non-equilibrium state of the particles, we will set some Cartesian 
coordinate systems at each point (the monad) inside the particle /. 


The coordinate system in each point which is rotation-fixed to the mass is 


recorded as Ko (Sats } , the coordinate system in the point (xx), whose 
coordinate axes are (nomm) . The coordinate system K! Ca) refers to 
the coordinate system in the point (%,,.x),.;) , the axis directions of which are the 
same as( 495.0 ) . The origins of both Ks) and K (3k) are the 
centroid ¢,(%,,x,,x;,/) of the mass of /(x,,x,,x,) in the point (the monad) 


(cides | 


Figure 10 Coordinate systems K,,K;,K, 


The axis directions of K! (aes are fixed in space, but it is in translation 
motion with the particle. However K; (Sata) is rotation-fixed with the mass of 
the particle / in the point (4535) , and rotates with the particle. The origin of the 


coordinate system K, (xxx) is (XX), a fixed coordinate system whose 
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axis directions are (x,,x,,x;) (See Figure 10). 

As shown in the figure 10, the angles between the axes of K; and K! are 
recorded as P, (xX Xbt',l ). It is the angle between the coordinate axis x; of 
K! and the coordinate axis 7, of the coordinate system K; ; 

Let COS P, = Q; (60) 
So one has A,A + AO + Ay, Ay, = Oy (k,l =1,2,3) (61) 
Therefore, only three of &, are independent. Give these three of &, and can get all 
nine Q. 

The coordinate system K; is rotation-fixedly connected to the material, so the 
rotation of K; is just the rotation of the mass. And let K; rotate round axis x, of 
the system K? , the rotation angle is ¢(%,,%,,%;,40',/) . 

Assume that 

Poy (x,.%/.%6051)-9,, (x; -Ep X Xbt) = Aq’, (otte) 

b, (xx xot t,1)-4,(x; eA ahta) = A¢' CEER 

(62) 

P, (dade al = Py (x,xX,x%,t,t'+dt', 1) =, (Seah!) 

Of, (XXX, t 1) = G (,%p,%5 00,0 + at, 1) — 6 (x,,%,,%5,0,0,0) 

(63) 

It is obviously that OP, can be obtained from (ò$ , OP, , õp ) . 


2, Force couple moment on the interface 


The elastic force action on the interface of two adjacent monads exerted on 


l (x1,X23x3) is reduced to one force (i.e. the aforementioned force 


F.(%,3%35X35t,t,/) ) plus one force couple moment, the component of which in the 
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axis x, is P(x. %extt']). 
Let the moment of force couple generated by the deformation of the particle / 
on the interface of two adjacent monads be proportional to Ag, . That is to say, 


there are 


Ti 


P Cee =—ah¢' (x... tt E 
E 


L 


=-a| ¢ (xox xt) -4 (x; — 614% )5%;st,t51) |— 
E 


L 


P, (X,.%).%pst.t',1)= add (x, + 6,X)%ptoll)— 
E 


i 


=al% (x; +E Xp tpb t I) -4 (xxt) 
E 


i 


=-a,| ¢ (ase aptt =o (a, E | 
ĉj 


=a,|¢ eae +e tph t )-4 (xx 601) -— 
j 


P- Tere =-a,Ag* (x,,2).%.t,1',/)— 
; E, 


=-a,| ¢ CERERE E Serpe 
Er 


Plett 1) akg Cee tapt tI) 
E, 


=a,|¢ (a em ee Le (64) 
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Here the subscripts to the couple moment p represent the surface (the former) and 


the component direction(the latter). Therefore B- oP. oP P. pPo pei are the 


components of x, -direction of the couple moment on the surfaces of 
i ,i*, j ,j ,k „k“ respectively . The right sides of the later four relationships in (64) 


are the curvatures of particle / (xx EA . The nature of (64) is the constitutive 


equation. 


3, About momentum moment theorem 
Write the momentum moment theorem for / Ce ree in the coordinate system 


K, (daa) i 


2 


M (x,t,t,1)& dt’ = L(x,t,t'+ dt',1)-L(x,t,,1) 


+E" (x,t,t > t'+dt',l)—L* (x,t,1' >t + dtl) 0?) 


Here M bethe force moment on / (xxx) , L bethe momentum moment of 


I(x,,%,,%3) 72" and g” bethe momentum moments of the outflow from and the 


xX 


Figure 11 Schematic diagram of momentum moment 


calculation 
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inflow into the monad (xxx) . These moments all have the reference point of 
(%,,%,,x;) , the origin of coordinate system K, (x,,x,,%;). 
In Figure 11, R is the vector of the centroid c, in the coordinate system 


K, , the unit vectors on three axes of the coordinate system K, are u, „Uz, lbs 
respectively. 


The momentum moment of l (yds) is 


2 


of a — o2 —_ 
L(x,t,t l)=e R(x,t,0,1)xU(x,t0,1)m(x.00,1)+ > € rx AmU! (66) 


o2 = SaN š r b 
Where, TP = mxr! , @=@(x,t,t',l) is the angular velocity of the particle 
U, =@xné 8 y p 


L(x X3x). 


By calculation there is 


2 


L) a R(x,t,0',1)xU (x,t,0/,0) m(x,t,0,1) 


+(1,@, —1,,@, —1,,0,)u, + (IQ, —1,@, -1,0 Ju, + (13, -1,,0, —I,,0,)u, (67) 


4 


4 
o oO 
12 12 Pagal 
Among them, /,, = > Am, (x2, +x0,)é » I= > AM, X; Xg E 
i i 


That are the moments of inertia of / (Rites ts) in the coordinate system K!. 


In the equation (67), @, = o,(x,t,t',1), Ly =I. (x,t,t',1) y damla (x,t,t', L). 


And let the masses flowing out (or into ) through the six interfaces be denoted as 


(4,4,4,4, As, A, ) respectively. The vectors of the centroids of these masses on the 


interfaces are CA (relative to the origin of the system K, ), or 


(hi, hs, hs, h, hsh). And that the h, is just (x A ae ), which have been already 


qi? q, j?” gk 


obtained in the previous section. Therefore these vectors are known or available. 


Then there is 
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a 6 o2 = 
L" (x,t, >t t+dt,l)= > R(x,t,0 De xdm" (x,t, >t +dt',1)U (x,t,t',l) 


q=l 
o4 


6 
DA x[ v (x6 t I) xri (x00) Je Sm (x,t,t! >t! +dt',1) 
q=l 


o2 
L" (x,t, >t +dt,1)= F R(x,t,t le xm” (x,t, >t +dt',1)U (x,t,t',1) 
q 


6 o4 
+17 (xt, 1)x] O(x, 0,01) x7 (x,t,t,1) |e Sm” (x,t,t' —>t'+dt',1) (68) 


q=l 
The three components of i (aay => t' + dt',1) are: Le" (x, t,t >t +dt',l) 


The three components of I” (x,t, ttt dt',1) are: LY (x,t,t' > t +dt',l) 


The three components of 5 RE t,t',1)xU(x,t,0', 1) m(x,t,t', 1) are L; (x,t,t',1) 


Then the momentum moment theorem can be written as 
ð 2 
M (x,t,t',l)e: dt’ =L (x,t,t'+ dt',1)- L? (x,t,0,0) 


+(112, -1 p0, -130 ) Ea -(1,2, -1 p0, -130 ) 
(t,t'+dt'1) 


(t,'./) 


+L" (x, t,t >t +dt',1)- L (x,t,t' > t'+dt',1) 


o2 
M, (x,t,t',l) a: dt' = L; (x,t,t'+dt',1)- L; (x,1,0,0) 


+(1,0, — 1,0, — 1,30, kaai -(I,0, -120 —1,,0,) 


(ttl) 


+L (x,t, t > ¢' +dt',1)-L} (x,t,t' > t'+dt',1) 


o2 
M, (x,t,t' l)e: dt’ = L} (x,t,t' +dt',1)- L (x,t,0',1) 


+ (1,30, — 1,0, — 1,0, Jetas - (130, -13o —1,,0,) 


(¢,t',/) 
+L (x,t, >t +d’ I-E (x,t, >t +de) m 
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Where the force moment M (x,t,1',/ ) is the force moment to be imposed on 
I (xX, X;) (relative to the origin of the coordinate system K, ), and 


M, (x,t,t',1) =P. (x,t,t',1)+ m}, (x,t,t',1)+m{,(x,t,t',1)+m,(x,t,t',1) (70) 


There are four items on the right side of the equation (70). The first item is the 
S — component of the force couple moment on the interface of the adjacent monads 


due to elastic deformation; The second is the $— component of the moment of the 


elastic force F, its action point is the centroid of the cross-section between the 
adjacent monads, referring to the origin of the coordinate system K,; The third is the 
s— component of the moment of the force f acting on the particle /(x,,x,,x;) 
due to the other particles in the monad (The action point of f is the centroid c, of 


/—mass in the monad.) referring to the origin point of the coordinate system K, ; 


The fourth term is the $— component of the couple moment acting on /— particle 


due to the other particles in the monad. 


4, Solving rotation problem 
Given ieee) alert eae i) 
By using successively the following steps, 


(1) One finds from ¢, that 
Ad: (cape 1) =¢, aaa —¢, (x; a) 
(2) By using of (64) one can getat P (x,t,1’,/) 
(3) Formula about the force moment is 
M, (x,t, [=P (x0) +m, (xttl) +m ,(x,0,0,1)+m, (x,2,01) (70) 
It can be seen from (69) that M, (x,t Aal ) ~é€ a high-order small quantity. 


Therefore, there is an approximate relationship. 


m, (x,t,0,1) =P, (x,t,t',1)-m} (x,t,t',1)-m{,(x,t,t',1) 
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Thus, take 


mi (x, 6,01) =—P (x,t,0,1)—m5, (x,t,t',1)-m{,(x,t,t',l) (71) 
Let m? be approximately m, and then correct m; to get the more accurate value 
of m,. 


(4) Amendment of m* 


When the interaction between the particles is a contact action, then the sum of the 


internal force moments of the system (refer to any point) is zero. Thus there is 


k a k 
EIR e T N EE m a Se (72) 


I= , A 

In square brackets, it is the moment of force Wiesel ) and the origin of the 
coordinate system K, is taken as the reference point. Yet the force Wile tai a! ) has 
been already found in the translation problem. 


Replacing m, in (72) bym* in (71), (72) will become 


ae 


SIRI R(x,tt I)e xf (x,t,t,1) Waal O Ge 


i= > A 


And to allocate ©, (x,1,2') based on the value of m;(x,1,t',/), that is, let 


k 
(a1) =) 9.07) (74) 
/=1 


k 
and m; (x,1,0',0) =D b., $2, = 
l=] 


Q, (x,t) 


Q, (x,t,t',1)= 5 


— A, (75) 


Here b, is a positive constant. 
Take m,(x,t,0',/) =m: (x,t,t',1)-Q,(x,40,0) (76) 


Obviously substituting the m,(x,t,2',/) in (76) into (73) there is 
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o2 = k 
Dd | R(x,t,t', De xf (60,0) | +9 m, (x,t,t',1)=0 (77) 


k 
=I s H 


(5) Substituting so obtained m, (x,t,t',1) into (70), the moment of force 


M,(x,t,t',/) can be calculated. 
And by using of the momentum moment theorem (69) one can find 
@,(x,t,t'+dt',1), 


(6) Finding by use of @, (x,t,t',/) 


2 


g (x,t,t' + dt’,1)= 9, (x,t,,)+ 2 dt'a, (x,t,t',1) (78) 
Og, = ¢, (x,t,t' + dt',1)-¢,(x,t,0,1) (79) 
Then you can find óP, (x,t,t',L) from ô$, (x,t, t',1) . The state of rotation of 


1(x,,%,%) is indicated by (5¢,,69,, ). 
There is @,,(x,t,t'+dt',/) =9,, (x,t,0,1)+ 09, (x.t,051) (80) 


(7) Starting from @, (x,¢,'+ dt’,1) ,% (x,t,¢'+ dt',1) ` Py (x,t,t'+dt’,/), repeat the 
above steps again. 
The moments of inertia of / (x,t,t') in the system K! are required in the 


calculation. This will be discussed below. 


5, About the moment of inertia 
The moment of inertia of the particle-mass m(x,t,t',/ ) in the coordinate 
system K; (x,t,t',]) ,recordedas J Ke (x,1,0',/), is given . The moment of inertia of 
m(x,t,t',) in the coordinate system K/(x,t,2',/), recorded as L (x,t,t', l), can 
be obtained by use of 9,,(x,t,t'/) and I9 (x,t,t',1). T; (x,t,ť,1) is just the 
moment of inertia in the momentum moment theorem (69). 
The expressions of 7 j (aia ) are given below. 


ol 


l o,l o,l 
I, = (ana, + atai + (ana, + apa Ie +(ana + RTA 
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-a a I? -a 


ol o,l 
a^ ji naala — 0 133 (81) 


Ty = Oly + Ogle Fada -atali Iita ~ tala (82) 
In the coordinate system K : (x,t, t'l ) let 
I (x,t,t' +dt' 1) =I! (x,t,t',1) 
+81 (x, t,t! >t +dt',1)- SI" (x,t, >t'+dt', I) (83) 
When the coordinate system K (x, t,t', d ) is translated from GC, (x,t „t'l ) to 
G; (x, t,t’ +dt',1), the coordinate system K; (x,¢,0',/) will be K; (x, t,t'+dt',1). After the 
coordinate transformation the value of / oe (x,t „t +dt',l ) will become I He (x, t,t'+dt',l ) in 


the coordinate system K; (x,t,t'+dt',l). Then using the angles P, (x,t,t'+dt',1) the 


moment of inertia I; (x,t,t'+dt', l) in the coordinate system K! (x,t,t' +dt',l) can be 
obtained. 

In solving the rotation problem, it is also assumed that the moment of inertia 
caused by the inflow and outflow masses is calculated in the case of the inflow and 
outflow masses being concentrated on there centroids on the interfaces. For example, 


the shaded portion in Fig.12. is the mass A which flows out of the monad 


ee, and enters into the monad (x; -Ep XX). The centroid of the mass on 


the interface is o(a!" =i a) in the coordinate system K (5a) , 


be ot X, 


l l L 
Figure 12 Finding the moment of inertia 


of the inflow (outflow) masses 
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and chx; sia a] in the coordinate system K, es ZEX Xpt, t1) . Due to the 
relationships between coordinate systems K, (lotat) , K! (aa AA and 
ER the coordinates of point € in K; are obtained. 


So by focusing the 4A mass on the centroid € , you can calculate the moments 


of inertia of 4A in the coordinate system HERTA and 
K! (x, -EX ee , which are OL (x, Xbt >t +dt',1) and 


ore" (a -En Xp Xp bht OE + dt',1) respectively. 


8 static state of particle system 
Static state means that all particles in the system are at rest. In this case, in order 


to understand the state of the particle system, it is necessary to understand the state of 
the forces of the particle system. The forces on /(x,t,¢") are: the elastic action on the 
corresponding interface (*-surface), the action can be divided into elastic force 


F.(x,t,t',1) (action point c'on x -surface) and elastic force couple moment 
P (x,t,t',1); the force to /(x,t,t') applied by other particles (or particle parts) in the 
monad (the force action point is c, ,the centroid of /(x,t,t')) f.(x,2,0',/) , as well 
as the force couple moment m, (x,t, t'l ); gravity F® (only the component in the 


vertical direction) , its action point is cĉ, too(see Fig. 13.). 


X—E X 


Figure 13 Forces on the particle /(x,t,¢’) at rest 
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Since /(x,t,t') is static, there are the balance equations 
F (x,t,t',1)+ f, (x,t,t',1)+ FE (x,4,0,1) =0 (84) 
(RxF(x,t,t',1)), AP (x61) +m, (xt) +R F (x01) HRF (x80) =0 (85) 
Here R,R’ are the vectors of c, , c' respectively. 


The referred point of the force moment is the origin of the coordinate system K,. 


k 
Because > f,(x,t1',/)=0, Sn x1,1)+ PERD) =0 (86) 


l= 


So, let the (84), (85) be summed over l(x,t,t"), it is obtained that 
VE (3,47, = 3 (ata) (87) 
I=1 il 
> (R R'xF(x,t,01)) + ÖP P (x,t,t',1)= -X (Rx FF wD (88) 
= i=l = 

The right sides of equations (87) and (88) are gravity and gravity moment 
respectively, while the left sides show the elastic force and elastic force moment 
generated on the interface between adjacent monads due to particle deformation. 
Therefore, equations (87) and (88) give the conditions that must be satisfied by the 
elastic force and moment generated on the interface between adjacent monads due to 


the deformation of the particle. On this basis, the static state of the particle system can 


be further explored. 


9 Conclusions 
(1) The basic idea is that the time-space in which there is the movement system of a 
large number of particles consists of the monads. The volume of a monad is infinitely 
small, but not zero. The physical quantities throughout the interior of the monad are 
not the same, i.e. are non-uniform. The physical quantities of the elements of the same 
particle in the monad are same, while that of different particles (or particle parts) in 
the monad are different, which are discontinuous. However, between the different 


parts of the same particle in the adjacent monads, the difference of the physical 
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quantities is infinitely small, which in turn reflects some continuity of granular motion. 


(2) The particle motion considered in one monad occurs within ~g¢? ofthe spatial 


range. This small scale range of motion gives a basic approximation in solving the 
problem. That is, the sum of the forces applied on the particles (or the particle parts) 
should be high order infinitesimal, the sum of the force moments applied on the 
particles (or the particle parts) should be high order infinitesimal also. 

(3) When the number of particles of small scale(« £) is large, for such a large 
number of small particles, we cannot use the method of this paper to solve one by one. 
One reasonable way to do is to find out the total effect of these many small particles. 
(4) The translational motion can be solved independently. After solving the 
translational problem the parameters on rotation motion can be obtained and the 


rotational problem will be solved easily. 


(5)The time step is taken as (t —t'+dt') in above discussing. People can divide 


whole time interval (t =0—t'=T) into ” segments. It is reasonable that the step 


istakenas At’ (i.e. nAt'=T ) in calculation. 
(6) On the initial value of the formulation, the following initial values should be given 
previously, 


for translational motion only 


U oner 01). Fae et =01), Mtaa st =h) 
(and its centroid position R(x, x,,x,t,t =0,1) ), the cross section s;,(x,4,t'=0,/) , 


S; (x-, t,t =0,1 ), the centroid-position of the cross section x a ae (t „t =0,1 ) 


if considering both translational and rotational motion, need an addition of initial 


values 


AEREE =0,/), Ole ei =O), Py (ae =0,/) 
bes ae =0,/) 


The conditions must be met between them that the differences of the different 


initial values of the same particle between adjacent monads are infinitely small. The 
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initial values can be selected under the condition that the compatibility is satisfied. 
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